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SUMMARY 

An analytical study was performed to provide information on the heat 
loss to space from three radiator finned-tube configurations in which 
there is radiant interaction between fins and tubes. Numerical results 
for the heat-transfer characteristics of these configurations were ob- 
tained under the condition that the surfaces are blackbody radiators. 

A comparison showed the heat-loss results for the various configurations 
to be of very similar magnitude for equal fin weight. A method of analy- 
sis for the case of nonblack diffuse surfaces was also formulated. 


INTRODUCTION 

The use of radiating fins as a means for augmenting the heat- 
transfer performance of space radiators has received wide consideration 
in recent years. In response to the need for design information, ana- 
lytical studies have been carried out to determine the rate of heat loss 
from fins to space under the combined action of surface radiation and 
internal heat conduction. The majority of such investigations has been 
concerned with fins that are essentially thin plates whose surfaces may 
radiate to space. Heat is transferred into such a fin by conduction 
through an edge that joins a tube within which flows a hot fluid. 

In the initial studies of such plate fins, a simple model.was 
adopted wherein the radiant interaction resulting from mutual irradi- 
ation between fin and tube was not included. References 1 and 2 are 
representative of several analytical treatments carried out along these 
lines. Radiators of practical interest, however, may consist of finned- 
tube configurations in which there is substantial radiant interaction 
between radiator elements. In practice, the detailed arrangement of 
fins and tubes in the radiator may be determined by considerations re- 
lating to fabrication, bonding, and assembly, to structural rigidity, 
and to meteoroid protection. Ultimate design decisions may be based on 
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a compromise among these several factors. It is therefore of interest 
to know the heat-transfer characteristics of a range of radiator con- 
figurations. 

The aim of this investigation is to provide basic heat-transfer 
information for three finned-tube configurations, all of which involve 
radiant interaction between elements of the radiator. The analysis has 
been carried out with the dual purpose of providing specific information 
for each of the configurations and also of showing comparisons among the 
heat-transfer performances of the several configurations. Numerical re- 
sults were obtained for the case in which the surfaces of the fins and 
tubes are blackbody radiators and the tubes are squared off. The method 
of analysis for nonblack surfaces is also formulated in detail. 

Studies of plate fin and tube configurations that include radiant 
interaction between radiator elements are reported in references 3, 4, 
and 5. Unfortunately, the results of references 3 and 4 have to be 
regarded with some degree of uncertainty. This is because the angle 
factors in the radiant interchange analysis were not correctly derived. 
In addition, these references formulate the nonblack case by simply 
inserting factors of emissivity e as multipliers of the emission 
and absorption terms that appear in the blackbody analysis. Such a 
treatment is only correct when there are no interreflections between 
the participating surfaces. When interreflections do occur, it is 
necessary to take account of the fact that the radiant flux leaving 
a surface location includes both reflected and emitted energy. Refer- 
ence 5 is concerned only with the central fin-tube configuration. 
Numerical results are given for the case of black radiating surfaces. 
These results are not precisely comparable with those of this study 
inasmuch as the present analysis employs a squared-off model of the 
tube surface. 

This research was sponsored by the National Aeronautics and Space 
Administration through the Office of Grants and Research Contracts. 


SYMBOLS 

A surface area 

B radiosity, radiation leaving surface per unit time and area 

F angle factor 

H incident radiation per unit time and area 
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"base surface height 
thermal conductivity of fin 
half-height of fin 

3 2 / 

conductance parameter based on dimension L, /kt 

rate of heat loss to space from fin and from portions of 
base surface that interact with fin 

absolute temperature 

temperature of base surface 

thickness of fin 

dimensionless coordinate, x/L 

position coordinate on fin 

dimensionless coordinate, y/L 

position coordinate on fin 

dimensionless coordinate, z/h 

position coordinate on base surface 

absorptivity 

dimensionless radiosity, B/crT^ 
emissivity 

temperature variable, T/Tb 
reflectivity 

Stefan-Boltzmann constant 
angles defined in fig. 2 


P 

a 

cp,cpW 

Subscripts: 

b base surface 

c due to conduction 
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f fin 

rad radiation 

s space 

1 base surface 1 

2 base surface 2 


GENERAL CONSIDERATIONS 

The specific finned-tube configurations that will be considered in 
this investigation are pictured schematically in figure 1. In the left 
diagram of figure l(a), the tubes are contained between a pair of plane, 
plate fins. This arrangement will be designated herein as the closed- 
sandwich configuration. In practice, it will neither be possible nor 
desirable to have a small region of contact between the tube and the 
fin. Such a limited contact would give rise to a large thermal resist- 
ance tending to inhibit heat flow from the tube into the fin. It would 
appear desirable to employ fillets of high-thermal- conductivity material 
to increase the area of contact. This suggests that a suitable model 
for analytical treatment of the radiant interchange within the space 
between the fins might be a squared-off tube. This idea has been adopted 
herein and is pictured in the right diagram in figure 1(a). The second 
configuration to be considered is similar to that already discussed, 
except that one plate is removed. This is pictured in figure 1(b), with 
the left diagram showing circular tubes and the right diagram the model 
adopted in the analysis. This assembly will be designated as the open 
sandwich. The last of the configurations that will be analyzed is the 
central fin and tube radiator, which is pictured at the left in 
figure 1(c). To permit meaningful comparisions with the prior cases, a 
squared-off model of the central fin will be employed as shown at the 
right in figure l(c). 

The length of the various finned-tube assemblies in the direction 
normal to the plane of figure 1 is sufficiently great so that end effects 
need not be included in the derivation of the radiation angle factors. 

temperature variations in this direction are also assumed 
sufficiently gradual so as to maintain the two-dimensional character of 
the problem. ~ 

Dimensional nomenclature and coordinates appropriate to the analy- 
sis are shown in figure 1. The squared-off tube surface, hereinafter 
called the base surface, is assumed to have a uniform temperature Tv,. 
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If consideration is given to the parameters that govern the heat trans- 
port, it is found that, after dimensionless variables have been formed, 

there emerge: 


L/h, (N c ) t = *tJl 2 


/kt 


( 1 ) 


The first of these is purely geometrical and governs the aspect ratio of 
the cavity formed by the fins and tubes. The N c parameter is essen- 
tially a ratio of thermal conductances, comparing the radiation capa- 
bility at the surface with the heat conducting capability of the m. 

For highly conducting fins, is very close to zero; for lesser con- 

ducting fins, N c takes on larger values. Hereinafter, N c will be 

termed the conductance parameter* 

The number of emissivities e and absorptivities a that enter 
the problem depends on the details of the particular system. The radi- 
ation properties of the base surface may either be different from or the 
same as those of the fins. There may or may not be a difference between 
the radiation characteristics of surfaces radiating to space and those 
radiating internally. If radiation from an external source such as the 
sun is involved, the absorptivity for solar radiation may be different 
from that for infrared radiation. The description of radiation from an 
external source would involve at least two additional parameters, one 
for the magnitude of the incoming energy and a second for its angle oi 

incidence. 

In view of the large number of parameters, it is clear that an 
attempt at a completely general treatment would involve a massive compu- 
tational effort, especially since numerical solutions of nonlinear, 
integrodifferential equations are required. This suggests that the role 
of each parameter be carefully reconsidered. In this connection, it is 
reasonable to expect that surface materials having the highest practi- 
cally obtainable emissivity will be used. Such emissivities may very 
likely be attained by the use of coatings. Thus, in practice, t e 
emissivity should not differ too much from that of a black surface. 

With this in mind, the main emphasis of this study will relate to blac 
surfaces, and this is the condition for which numerical results will e 
obtained. The use of black surfaces results in a computational problem 
much more tractable than that encountered when nonblack surfaces are 
used. However, for completeness, a formulation of the problem for non- 
black surfaces will be carried out in the appendix. In addition to the 
condition of black surfaces, the present analysis will also be confine 
to the situation in which radiation from an external source is negli- 
gible compared with energy emitted by the finned-tube system. 


6 


TEMPERATURE DISTRIBUTION 
Governing Equation 

In carrying out the analysis. It will be convenient to begin with 
the closed- sandwich configuration shown in figure l(a). Once the 
governing equation for the temperature distribution has been formulated 
for this case, only simple modifications need be made to obtain the cor- 
responding equations for the other configurations* 

Close d -sandwich configuration . - The analysis begins with the appli- 
cation of the conservation-of -energy principle to an elementary volume 
of fin, which is shown as shaded in figure 2. Under steady-state con- 
ditions, the net outflow of energy from such an element must be zero: 
therefore, 

(dC! net> c + (d W ra4 = 0 (2) 

For a thin (one-dimensional) fin, the net outflow of heat due to con- 
duct ion is 


(dQ net> - ^ foe (3) 

G dx 

per unit length normal to the plane of the figure* The net outflow of 
radiant energy is the difference between the emitted and the incident 
radiation* The radiant emission is given by the Stefan-Boltzmann law as 

[crT (x)J 2 dA x , in which the factor 2 is introduced because both surfaces 
of the element are emitting. With this, the net radiation may be 
written 

( dQ net) rad = [2 cjT 4 (x) - H(x)] d^ (4) 

where H denotes the incident energy falling upon x per unit time 
and area* 

The incident radiation comes from two distinct zones: (l) from the 

opposite fin, and (2) from the base surfaces* The contribution from the 
opposite fin will be derived first. Consider radiation leaving a typi- 
cal surface element of the opposite fin, say dAy* This element is in 
reality a narrow strip, having a dimension dy in the plane of 
figure 2, but extending to infinity in both directions normal to the 
plane of the figure. The energy emitted by dA y is ^aT 4 (y)JdA y , 

and this is distributed over all directions in space. Of this, an 
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amount crT 4 ( y ) cLAy dF 


is incident on dA^., where dF^..^ 
an angle factor for diffuse interchange. However, by the reciprocity 

0 -t ft -I TH T.T-? 4- V, 


dAy dA x 


IS 


rule for diffuse angle factors, cLA^ dF^ -dA = dA^. 

J V x 


dF 


dA -dA * 
x y 


this, the contribution of 
-4 


With 
is 


y x 

cLAy to the incident energy at dA x 
crT 4 (y)dA x < ^’dA x dAy^ ^ut such contributions arrive from all locations on 
the opposite fin, and the total amount is found by integration: 


/ ^ 4 (y)dF a , „ (5) 

f y 

Jy=0 

The radiation incident at dA x from the base surfaces may be de- 
rived in a similar manner. The final expressions are most easily ob- 
tained by modifying equation (5). Since the base surface temperature 
is uniform, it may be removed from under the integral sign. Additionally, 
the integral of an infinitesimal angle factor yields a finite angle 
factor. Applying this to the two "base surfaces 1 and 2 (fig. 2) gives 

"x^( F 4A x -1 + ? dA x -2) (6) 

With expressions (5) and (6), the radiant energy H(x) incident per 
unit area at x is 


/ 2L 

+ oT b( F dvi + X- 2 ) (? 

. 0 

The net radiation may then be found by substituting H(x) into equa- 
tion (4). Thus, with both the net radiation and net conduction avail- 
able, the energy balance equation (2) may be evaluated. This gives 

/ 2L 

oT 4 (y)dF dA _ dA - ^b( F dA -1 + F dA - 2 ) 
x y V x x / 

0 


( 8 ) 
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Introducing the dimensionless variables 

X = x/L, Y = y/L, 9 = T/T t (9) 

equation (8) becomes 



( 10 ) 

in which the conductance parameter N c has already been defined by 
equation (l). Because of the symmetry of the problem, it is clear that 
®(X) and 9(Y) are the same functions; only the independent variables 
have been interchanged. Thus, equation (10) is an integrodiff erential 
equation for 0. Since both 9 and 9 ^ appear, the equation is non- 
linear. 

To complete the specification of the problem, the boundary con- 
ditions must be stated and expressions must be provided for the angle 
factors that appear in equation (10). One of the boundary conditions is 
that T = T-^ at x = 0. Additionally, there is a symmetry condition 
at x = L, the mathematical representation of which is dT/dx = 0. In 
dimensionless terms, these are 

0(o) = i (d 0 /dx ) x=1 = o (11) 


In the derivation of the angle factors. It is advantageous to make 
use of special properties associated with surfaces that have one elon- 
gated dimension (essentially infinite normal to the plane of fig. 2). 

In general, the angle factor between an element dAy and an infinite 
strip dAy is (ref. 6) 

^dA^-dAy = £ < sin ^ (12) 

wherein the angle cp is defined by the normal to dA and the con- 
necting line between dA^. and dA as shown in figure 2. From the 
geometry of the figure 

sin cp = (y - x)/\/(y - x) 2 + h 2 



9 


Carrying out the indicated differentiation results in 




For interchange between cLA^ and an infinite strip that spans a finite 
dimension in the plane of the figure, a generalization of equation (12) 
may be applied. For instance, for the angle factor to the base surface 
2 , 

F^ _ 2 = - ( sin 9" " sin 9' ) ( 14 ) 
x 2 


(Y - X) : 


+ 


h\2 

L. 


3/2 


— dY 


( 13 ) 


where cp" 
equation 


and cp' are defined in the lower sketch of figure 2. When 
(14) is evaluated and dimensionless variables are introduced, 


C ft <- 


7 


i - 


( 2 - *). 


(2 - X) 2 + (I)' 


q/2 f 


(15) 


Similarly, for the base surface 

1, 

r 



3 


X 


F dA x -l - ?< 


’** + (# 

T 72 r 


L 


J 


(16) 


With the foregoing formulation of the angle factors, the evaluation 
of the energy conservation principle is now complete. There is, however, 
one additional simplification that may be introduced into equation (10 ). 
This results from the symmetry of the temperature distribution about 
x = L; that is, T(x) = T(2L - x). When this symmetry property is in- 
corporated into equation (10) along with the appropriate angle factors, 
there follows 



r 

X 


' [ ^ (l/ 

VT 

1/2 




2 - X 


’< 2 - X) 2 * (if 

1/2 


J 


( 1 ?) 
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Inspection of equation (17) reveals that there are two independent 
parameters that must be specified for each solution; the conductance 
parameter N c and the aspect ratio h/L. Further study indicates that 
analytical solutions cannot be achieved, and numerical means must be 
employed. The numerical solutions will be discussed after the governing 
equations for the other configurations have been derived. 

Open-sandwich configuration . - The energy equation for the open- 
sandwich configuration is most conveniently derived by modifying the 
foregoing analysis for the closed sandwich. In fact, for the open sand- 
wich, it is only necessary to delete the terms that describe the radiant 
interaction between the fins - the integral terms appearing in equation 
(17). With this modification, the governing equation for the fin- 
temperature distribution in the open-sandwich configuration becomes 




20 4 (X) - 1 

r 





The boundary conditions remain as stated in equation (ll). 

By comparing equation (18) for the open sandwich with equation (17) 
for the closed sandwich, a very significant mathematical difference is 
apparent; namely, the former is a differential equation, while the lat- 
ter is an integrodifferential equation. In general, differential equa- 
tions are easier to solve than integrodifferential equations. While 
equation (18) is, in fact, simpler than equation (17), numerical solu- 
tions are still required. These will be discussed later. 

Central-fin configuration . - The energy equation for the central 
fin may, in turn, be derived by modifying equation (18) for the open 
sandwich. The aspect requiring modification is the interchange between 
the fin and base surfaces. For the central-fin configuration, the inter- 
change occurs between dA x and a wall of height h/2 (assuming t < < h) 
rather than of height h as in the open sandwich. Thus, wherever h 
appears in equation (18), it is replaced by h/2. Additionally, both 
the upper and lower sides of the central fin exchange heat with the base 
surfaces, which requires that a multiplying factor of 2 be inserted in 



the angle factor terms of equation (18). When these modifications are 
made, the revised form of equation (18) now applicable to the central- 
fin configuration is 



The boundary conditions of equation (ll) continue to apply. Numerical 
solutions of equation (19) will be discussed in the following paragraphs. 


Solutions 

Method. - As has already been indicated, it was necessary to use 
numerical techniques to solve the foregoing equations for 9 . The nu- 
merical solutions were carried out on a Univac 1103 electronic digital 
computer. The procedure for solving equation (17) will be described 
first. In essence, a weighted iteration technique was employed. The 
first step was to prescribe the parameters N c and L/h; then, a trial 
distribution of 0(Y) was selected. For a fixed value of X, the inte- 
gral appearing on the right side of equation (17) was carried out 
numerically. Then, another value of X was chosen and the integral 
was carried out once again, and so forth. In this way, corresponding 
to the particular 9 function, a value of the integral was generated 
at every mesh point X in the range 0 < X < 1. This set of values 
may be denoted by f (X), while the remaining terms on the right side, 
aside from 0 4 , maybe denoted by g(x). Thus, equation (17) maybe 
written as 


= (n_) [2 e 4 + f(x) + g(x) 

dx 2 L L 


( 20 ) 
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Integration of equation (20) subject to the boundary conditions (ll) 
leads to 


e(x) = 1 + 2 (n ) 

^ L 





e 4 



dX 



As previously described, f(x) is available for a particular trial distri- 
bution of 9 , while g(x) is a known function of X* Therefore, the 
integration indicated in the last term of equation (21) can be carried 
out numerically. Also, by utilizing the trial distribution of 9 } the 
integrals involving 9 may also be carried out. By these calculations, 
a new 9 distribution is generated from equation (21). 

In a direct iteration scheme, the new trial function would serve as 

input to the next cycle of calculations. However, in a nonlinear system, 

direct iteration may often lead to divergence rather than convergence. 

The procedure used herein was to take the input and output 9 functions, 

which may be respectively denoted by 9 and 9 and form their dif- 

^ n n*tl 

ference. The trial function 0 for the next cycle of calculations was 
then constructed as 

0* = 0 n + MVl " 0 n ) 0 < X < 1 


The size of A could be controlled manually from the console of the com- 

r 1 

puter. An on-line typewriter printed out / (0^ - 0 n )dX, and adjust- 

J 0 


ments could be made in A if the trend of this integral indicated di- 
vergence. In general, cases characterized by small L/h could be run 
at A values near 1, while cases at large L/h required smaller A 
values to ensure convergence (as low as A = O.l). 


Solutions were obtained for L/h values of l/2, 1, 2, 3, 4, and 5 
for each of which parametric values of N c were selected in the range 
from 0 to 10. In all, 56 cases were solved for the closed— sandwich con- 
figuration. 



13 


The same basic procedure was employed to find solutions of equa- 
tions (18) and (19), respectively, for the open sandwich and the central 
fin. However, there was one major simplification that resulted from the 
fact that these equations are differential rather than integrodif fer- 
ential; namely, the f(X) no longer appears in equation (21). The fact 
that it was no longer necessary to compute f(x) enabled the solutions 
for equations (18) and (19) to be carried out more rapidly than those 
for equation (17). The same parameter ranges were covered in the solu- 
tions of equations (18) and (19), as have already been discussed in con- 
nection with the solutions of equation (17). 


Results. - Each solution of the energy equations (17) to^(l9) P r °~_ 
vides a distribution of temperature along the fin. In view of the large 
number of cases that has been considered, it is neither possible nor 
necessary to present temperature profiles for all. Instead, a repre- 
sentative set of profiles will be presented for each of the radiator 
configurations. These results are plotted in figure 3 for the close 
sandwich, open- sandwich, and central-fin configurations. On each figure 
curves are given for T/% as a function of position x/L for para- 
metric values of and L/h. 


Figure 3(a) shows that the temperature drop along the fin is very 
small when N c is very small (see, e. g., the curves for = °: 25 J* 

This is due to the low thermal resistance of the fin. Additionally, . 
there is very little effect of L/h when N c is small, although it is 
clear that the temperature drop is greater as L/h increases. As c 
increases, so also does the temperature drop along the fin. The effect 
of L/h also becomes much more marked at the higher H c values. For 
a given value of N c , it has been verified that the temperature profile 
for the single, noninteracting fin falls slightly below that for 
L/h = 5. This is fully plausible, since the noninteracting fin is 
equivalent to the case of L/h = °°. 


The trends just enumerated for the closed sandwich are also in evi- 
dence in figures 3(b) and (c), which correspond respectively to the open 
sandwich and the central fin. The curves on these latter two figures 
fall somewhat lower than those of figure 3(a) because, for a given N c 
value, there is twice as much fin volume in the closed-sandwich con- 
figuration as there is in the open sandwich or the central fin. 


The figures generally show that, when N c is high, the tempera- 
tures in the central region of the fin are quite low. When it is re- 
called that the emission is proportional to the fourth power of the em 
perature, it is clear that the fin material in the central region is of 
marginal utility. The presence of a larger base surface (smaller L/h) 
tends to raise the temperature in the central region, making this par 
of the fin more useful. 
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HEAT TRANSFER 

Heat-t ransfer results based on the foregoing numerical solutions 
will now be presented and discussed for the separate configurations* 
Comparisons among these results will then be made* A superposition 
calculation of the heat loss will also be explored* 


Heat-Transfer Characteristics 

The heat-transfer results to be presented here will include the 
heat loss to space from the fin and from those portions of the base sur- 
face that interact radiatively with the fin* The heat loss from those 
po^lions of the base surface that radiate directly to space without 
interaction will not be included* In this latter category fall the upper 
and lower surfaces of the squared-off tubes, that is, the horizontal sur- 
faces of the squared-off tubes as they are pictured in figure 1* The 
heat loss from such surfaces is simply crT^ times their surface area, 
and this may be added at any time to the results presented herein* 


configuration . - When excluding the noninteracting 
portions of the squared-off tube as discussed previously, the heat loss 
to space from the closed-sandwich configuration is due to radiation from 
the surface of the fin. The rate of heat loss from the fin surface per 
unit area at x is cxT 4 (x). The overall rate of heat loss Q from the 
pair of fins that forms the sandwich is 

/ 2L 

oT 4 (x)dx (22) 

per unit length normal to the plane of figure 1. In dimensionless terms, 
this becomes 


— ^-r = / 0 4 (X)dX (23) 

««*b / 

•'0 

in which the symmetry of 8 about x = L has been taken into account. 

The group iLcrT-^ is the heat loss from the fins under the condition 
that their temperature is uniform at the value T^. Such a condition 
can be achieved only by fins of infinite thermal conductivity (N c = 0). 
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Therefore, 4LoTv represents the largest possible heat loss from the 
system. The ratio expressed by equation (23) thus compares the actual 
heat loss with the ideal heat loss. Such a ratio may be regarded as an 
effectiveness of the system. 

The effectiveness ratio represented by equation (23) has been 
evaluated by numerical integration of the previously discussed tempera- 
ture solutions. The results thus obtained are presented in figure 4 as 
a function of N c for parametric values of the aspect ratio L/h. In- 
spection of the figure reveals that for any fixed L/h the heat loss, 
from the system decreases with increasing values of N c . This. trend is 
in complete accord with physical reasoning, since an increase in N c is 
an indication of an increase in the thermal resistance of the fin 
(because of a decrease in thermal conductivity k or fin thickness t). 
The sharpest decrease in fin heat loss occurs at small values of N c . 

It may additionally be seen from the figure that the effectiveness 
decreases monotonically with increasing L/h. Since L appears in. the 
ordinate variable as well as in L/h, this trend may be most convenient y 
discussed by supposing L to be fixed. Then, an increase. in L/h cor- 
responds to a decrease in the separation h between the fins. In par 
ticular, the case of L/h = «. corresponds to the situation in which the 
separation h vanishes, and the pair of fins forming the sandwich co- 
alesces into a single fin of thickness 2t. Thus, the heat-loss results 
for the L/h = 00 case are identical to those for a single, isolated fin 
of thickness 2t. These isolated fin results were taken from reference 5. 

To explain why the heat loss decreases with decreasing ^h (fixed L), 
it is first necessary to describe the role of radiation within the 
cavity of the closed-sandwich configuration. This radiation provides a 
means of transport between the fins and the base surfaces and between 
high- and low-temperature portions of the fin, which is in addition to 
the heat conduction in the fin. In effect, the internal radiation and 
the conduction are parallel heat-flow paths. Thus, the presence of the 
internal radiation is equivalent to an increase in the thermal conduc- 
tivity of the fin. The most important transport of radiation within the 
cavity takes place between the base surfaces and the fins. Either by 
physical reasoning or by inspection of equations (15) and (16), it is 
easy to see that the transport between fin and base surface decreases 
as the spacing h decreases. Additionally, from equation (13), it is 
also seen that the transport between widely separated locations on the 
fins (i.e., Y substantially different from X) decreases as h de- 
creases. Thus, internal radiation becomes less effective as a transport 
mechanism as h decreases. The trend among . the curves of figure 4(a) 
showing Q decreasing with decreasing h (fixed L) is thus made 
plausible. 
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It is interesting to inquire about the quantitative relation be- 
tween internal radiation and fin heat conduction, the sum of which 
constitutes the total heat loss. The heat conducted into the fin at the 
base Q c was calculated by Fourier's law: 

Q c = -kt(dT/dx) x=0 (24) 

This equation was placed in dimensionless form and then evaluated from 
the solutions of the energy equation (17). The ratio of Q c to the 
total heat loss Q has been constructed, and the results are plotted 
in figure 5. It is clearly seen from the figure that the importance of 
heat conduction in the fin relative to internal radiation decreases as 
L/h^ decreases. This is in accordance with the discussion of the pre- 
ceding paragraph. Additionally, it is seen that heat conduction is the 
dominant heat -transfer path at low values of N c and that internal 
radiation grows increasingly important as N c increases. This follows 
from the fact that small values of I\I C indicate small thermal resist- 
ance for heat conduction in the fin. 

Open -s andwich configuration . - From a consideration of the right 
diagram of figure 1(b), it is clear that there is a heat loss to space 
both from the fin and from the vertical base surface walls that inter- 
act with the fin. The net heat loss from the fin (both sides) is equal 
to the heat conducted into the fin. From Fourier's law 

Q f = 2[-kt(dT/dx) x=0 ] (25a) 

in which the factor of 2 is included because there is an inflow at both 
x = 0 and at x = 2L. In dimensionless terms, this is 


Q f = 1 ( d£\ 

4LoT^ 2 (N C ) \ ^/x=0 
D L 


(25b) 


The net heat loss from the base surface is equal to the energy 
emitted by the surface minus the incident radiant flux. The emission 

is simply crT-^ h . per unit length normal to the plane of the figure. The 
incident radiation originates at the fin and at the opposite base sur- 
face. In mathematical terms, the net base surface heat loss is 
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in which the factor of 2 is included, because two base surfaces are in- 
volved. By introducing dimensionless variables, this becomes 


% 


z & - h- 2 : 



0 4 (X)F dAx _ 1 dX 


(26b) 


The angle factor F aA^_i is already available in equation (16), while 
the angle factor F^_2 for interchange between two parallel surfaces 
follows from reference 6 (eq. (31-49)) as 


F 


1-2 



2L 

h 


(27) 


The total heat loss Q from the system is found by adding and 

Q f : 


Q = ^ (28) 

4LoT 4 4LaT 4 4LcrT 4 

As written in equation (28), Q is expressed as a ratio with the 
quantity 4LoT 4 as with the closed sandwich. It will now be shown that 

4LoT 4 is precisely equal to the maximum possible heat loss from the 
system. For the case of infinite thermal conductivity (N c = 0), the fin 
temperature is uniform; that is, 0=1 throughout. In this instance, 
q can be evaluated from equation (26b) by removing 0 from under the 

integral sign. Additionally, 
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where is the angle factor from the fin to the base surface 1. 

Thus, 


% 

4 

4Ix7TJ 

D 


= — (1 
2L V 


f i- 2 ) 


- F 


f-1 


( 29 ) 


But, hy angle factor algebra, F 1-2 + F 1 _ f + F-^g = 1, in which F ] __ s 
indicates the angle factor from the base surface 1 to space. Clearly, 
because of symmetry, F 1-s = Fp.f. Also, from the reciprocity theorem 

for diffuse angle factors, 2LFj._^ = hFj_f. Incorporating this into 
equation (29), there follows 


4LcT^ 

D 


h_ 

2L 


F 


1-f 


(30) 


It is to be remembered that this applies for the limiting situation 
where the fin is isothermal (N c = 0). The fin heat loss for this 
limiting situation cannot be calculated from equation (25a) because k 
is infinite and d0/dX is zero* Instead, is determined as the dif- 

ference between the emitted and incident radiation. For the fin at uni- 
form temperature this is 

Of = 2 ^2LaT^ - hoTjp i _ f j (31) 

or 


Qf 

4IrfP< 


1 


_ IT 

2L 1-f 


(32) 


Now, adding equations (30) and (32) to find the total rate of heat loss 
gives 


Q = 4LctT£ (33) 

This is the maximum possible heat loss from the system. Thus, the ratio 
expressed by equation (28) is the effectiveness of the system as in the 
case of the closed sandwich. 

Values of the effectiveness have been determined by numerically 
evaluating equations (25b) and (26b) in conjunction with the 9 solu- 
tions of equation (18). The results thus obtained have been plotted in 
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figure 6 as a function of N c for parametric values of the aspect ratio 
L/h. The general trends that have been discussed in connection with 
figure 4 for the closed-sandwich configuration are also in evidence in 
figure 6 for the open-sandwich configuration. However, there are certain 
interesting differences in detail. The decrease in Q with increasing 
L/h (decreasing h) is much sharper for the open sandwich than for the 
closed sandwich. In considering this result, it should be noted that 
for a given fin thickness t there is only half as much fin volume in 
the open-sandwich configuration as there is in the closed-sandwich con- 
figuration. Consequently, there is less cross-sectional area for heat 
conduction in the open sandwich. It is therefore reasonable that the . 
open sandwich should be more sensitive to the effects of internal radi- 
ation than the closed sandwich. Consequently, as the internal radiation 
diminishes with decreasing h, the heat loss Q will be more affected 
for the open sandwich than for the closed sandwich. Comparison of 
numerical values between figures 6 and 4 will be reserved until later. 


It is also of interest to know what fraction of the overall heat 
loss is transferred from the fin and what fraction is from the inter- 
acting portion of the base surface. This information is obtained di- 
rectly from equations (25b) and (26b) for N c > 0; for N c = 0, equa- 
tions (30) and (32) are used. A plot of Qf/Q is presented in figure 7 
Corresponding results for Q b /Q may be readily obtained by applying the 
relation Q b /Q = 1 - (Qf/Q). Inspection of the figure reveals that the 

fin's fraction of the total heat loss is greatest for large values of 
L/h and small values of N c . The trend of Qf/Q increasing with in- 
creasing L/h is related to the increase in fin surface area relative 
to the base surface area. On the other hand, the changes in Qf/Q with 
N c are due to changes in fin temperature distribution. For small 
values of N c , temperatures throughout the fin are relatively high, and 
the fin heat loss is correspondingly high. For high values of N c , the 
fin temperature drops off rapidly with increasing distance from the base 
surface, and, as a consequence, the fin heat loss is low. 

The open-sandwich configuration pictured in figure l(b) is unsym- 
metrical in that the lower surface is a plane sheet, while the upper 
surface is, in essence, a cavity. The vertical walls of the cavity are 
at the temperature T b , which is the highest temperature in the system. 
It would therefore be expected that the heat loss to space from the 
cavity side would exceed that from the plane side. It is of interest 
to explore the relation between these heat losses quantitatively. The 
heat loss to space from the plane side of the configuration is given by 


^plane side 



crT 4 (x)dx 


(34) 
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Numerical results have been obtained by evaluating the integral with the 
9 solutions of equation (18). The ratio of the heat loss from the plane 
side to the total heat loss Q has been constructed, and the results 
are presented in figure 8 as a function of N c for parametric values of 
the aspect ratio L/h. As expected, the heat loss from the plane side 
is less than one-half of the total (except in the limit at N c = 0). 

From the figure, it is seen that, as the L/h ratio increases, the heat 
losses from the plane and cavity sides tend to equalize (i.e., the ratio 
plotted in fig. 8 approaches more closely to one-half). This is due to 
the fact that the base surfaces play a less important role as L/h in- 
creases (h decreases). Further study reveals that the unbalance in 
heat loss between the cavity and plane sides is accentuated at high 
values of N c . Large values of N c are effective in diminishing the 
fin heat loss relative to the base surface heat loss. Thus, at large 
N c , the cavity side is favored by virtue of the presence of the base 
surfaces. 

Central-fin configuration . - In common with the open sandwich, the 
heat loss to space from the central-fin configuration takes place in 
part from the fin and in part from the base surfaces. The net heat loss 
from the fin is expressed as before by equation (25b). On the other 
hand, the net base surface heat loss as given by equation (26a) must be 
somewhat modified to take account of geometrical differences between the 
open-sandwich and central-fin configurations. Figures 1(b) and (c) show 
that the radiant interchange in the central-fin configuration involves 
base surface elements of length h/2 instead of length h as for the 
open sandwich. Additionally, there are now four participating base sur- 
face elements rather than two as in the previous case. When these modi- 
fications are incorporated into equation (26a), it is found that the 
dimensionless base surface heat loss is expressed by 
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The overall rate of heat loss from the system is found by adding 
the contributions of the fin and of the base surfaces. Numerical values 
have been obtained by evaluating equations (25b) and (35a) in conjunc- 
tion with the temperature solutions of equation (19). The heat-loss 
results thus calculated have been plotted in figure 9 as a function of 
N 0 for parametric values of the aspect ratio L/h. The overall 
appearance of this figure and the general trends among the results are 
quite similar to those already enumerated and discussed in connection 
with figure 6 and, therefore, need not be repeated. A numerical com- 
parison of the heat-loss results of figure 9 with those of figures 6 
and 4 will be made in the forthcoming section. 

Inasmuch as there is a heat loss to space from both fin and base, 
it is of interest to inquire about the relative importance of each. The 
fin heat loss has been calculated from equation (25b), and, with this 
information, the ratio Q f /Q has been constructed. The results are 
plotted in figure 10, where it is seen that the fin's contribution to 
the overall heat loss is greater for large L/h and for small N c . The 
explanation of these trends has already been discussed in connection 
with figure 7. A numerical comparison of figures 10 and 7 shows that, 
for large L/h, there is little difference between the Q f /Q values for 
the open sandwich and for the central fin. This is fully plausible mas 
much as the base surface plays only a minor role when L/h is large. 

The detailed arrangement of fin and base surface should therefore have 
little effect. On the other hand, there are larger departures between 
the two sets of results at small L/h, with the Q f / Q values for the 
open sandwich (fig. 7) being higher than those for the central fin 
(fig. 10). 


Comparison on Basis of* Equal Fin Weight 

The foregoing presentation of results has been focused on the heat- 
transfer characteristics of each of the specific configurations under 
study. Now, consideration will be given to a comparison of the heat- 
transfer performance of these configurations. As a basis for the com- 
parison, the condition of equal fin weight has been selected here. 

Other bases of comparison for the complete radiator configuration are 
possible, one of which might be equal radiator (tube plus fin) weight, 
but this involves variables such as tube wall thickness that are not 
among the parameters of this study. 

Referring to figure 1, it is clear that the condition of equal fin 
weight is achieved when the thickness of eacn of the fins comprising 
the closed sandwich is one-half that of the fins in the other configu- 
rations. The fin thickness enters the analysis only in the^ N c parame- 
ters (see eq. (l)). Thus, for a comparison at equal fin weight, the 
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parameter for the closed-sandwich configuration should he taken 
twice as large as the N c parameter for the other two conf igurations. 

A comparison of the heat-loss characteristics of the various fin- 
tube configurations is presented in figure 11. This figure is sub- 
divided into several graphs, with each graph corresponding to a fixed 
fin weight. There are three curves in each graph to represent the re- 
sults for each of the three coni igurat ions. Figure 11(a) contains com- 
parisons for the lower N c values, while figure 11 (b) is for the higher 
N c values. The abscissa variable for all the graphs is L/h. 

Inspection of the figures shows that the major differences in heat- 
transfer performance among the various configurations occur at small 
values of L/h. This occurrence may be understood by recalling that the 
importance of the base surface heat loss is greatest when L/h is small. 
Therefore, it is natural that the details of the base surface arrange- 
ment should be important for small L/h. As L/h increases, it is seen 
from the graphs that the spread among the curves diminishes. In practice, 
it would be expected that L/h values in excess of three would be most 
probable. In such an event, the performance differences among the con- 
figurations are quite moderate* For example, in the middle graph of 
figure 11(a), which represents a practically realistic level of N c , the 
overall spread among the curves is about 3 percent at L/h =3. 

Generally speaking, the spread among the curves grows larger on a 
percentage basis as the level of N c increases. This is related to the 
increasing role of the base surface with increasing N . The N values 
°-^ figure 11(a) lie in the range of current practical interest, while 
those of figure 11(b) are somewhat on the high side. 

The open-sandwich conf iguration appears to show the best heat- 
transfer characteristics among the various configurations studied here. 
However, except at low values of L/h, the margin of difference is not 
decisive. 


Superposition of Results for Noninteracting Surfaces 

The heat-transfer results that have been presented in this report 
have been derived from rather lengthy numerical solutions. It is worth- 
while inquiring whether a close approximation to these results can be 
achieved by a combination of other, more easily derivable results. In 
this connection, it is natural to take as a starting point the heat-loss 
results for the single, noninteracting fin. These results are readily 
available in the open literature (e.g., in refs. 1 and 2). To these 
results for the noninteracting fin should be added a quantity that char- 
acterizes the heat loss of the base surface and, In particular, that 
part of the base surface that interacts with the fin. (As previously 
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noted, those portions of the squared-off tube that do not interact with 
the fin have not been included in this study. ) Perhaps the most reason- 
able heat loss to take for this purpose is that portion of the base- 
surface emission that streams directly out to space. For example, if 
the base surface emission for the open-sandwich configuration (fig. 1( )) 

is 2 ^oT-^h), then the portion that would go directly to space is 



where Fq__ s is the angle factor from one of the base surfaces (say, 
base surface l) to space. 

In the closed-sandwich configuration (fig. 1(a)) those parts of the 
base surface that interact with the fins are internal to the system. 

So, for this case, F 1-s would be identically zero. Thus, the super- 
position discussed previously cannot be applied to the closed-sandwich 
configuration. 

For the open-sandwich configuration, the base surface contribution 
Q* to the superposition is given by equation (36). In dimensionless 
form, this becomes 

Q > Lt < ■ (£) F 1 -s 

By angle factor algebra, F]__ s is related to the 
equation (27), as follows: 

*l-s *i- 2 ) 

Equation (36a) is to be added to the results for 
acting fin that are represented by the L/h = °° 

The actual superposition has been carried out in 
results are shown only for selected values of L/h to preserve clarity. 
The solid curves are the results that previously appeared in figure 6, 
while the dashed curves are those obtained by the superposition. The 
figure indicates that the heat-loss results predicted by the super- 
position method are, on the one hand, low at large values of N c and 
are, on the other hand, high at small values of N c - Indeed, for 
N = 0, the superposition method would predict an ordinate value that 
would exceed unity by an amount equal to the right side of equation (36a) 
Further, there is no reason to believe that, for L/h values exceeding 
5, the results from the superposition method may not fall high even for 
large N . From an overall appraisal of figure 12(a), it may be 


Fq_g, as given by 


(37) 

the single, noninter- 
curve of figure 6. 
figure 12(a), where 
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concluded that the superposition method applied here is not satisfactory 
for practical calculation purposes, at least for the open-sandwich con- 
figuration. 

The . same type of superposition has been attempted for the central- 
fin configuration. The base surface contribution is given by equations 
(36a) and (37) as before, but now the appropriate expression for 
F l-2 is equation (35b). A comparison of the heat-transfer predictions 
from the superposition method with those of the numerical solutions is 
shown in figure 12(b). The general characteristics of the comparison 
are very similar to those of figure 12(a), except that now the dashed 
curve lies above the solid curve for L/h = 5. The general conclusions 
relating to the practical utility of the superposition method remain as 
previously stated. 


CONCLUDING REMARKS 

The foregoing analysis has provided results for radiator finned- 
tube configurations in which the tubes have been idealized as having 
squared-off exterior surfaces. It would be expected, however, that the 
squaring off will not have a decisive effect on the heat -transfer re- 
sults in the L/h range of practical interest (i.e., larger L/h). 

An important contribution of this study relates to the conclusions 
that follow from comparisons among the results. Such comparisons on the 
basis of equal fin weight show that, for aspect ratios and N c values 
of practical interest, the heat-transfer performance does not differ 
appreciably among the various configurations that have been studied. In 
the practical range of L/h, it is not expected that the essence of this 
conclusion will be altered when round tubes (rather than squared-off 
tubes) and highly emissive nonblack surfaces are used. Under the con- 
dition of equal tube weight, it follows that optimizations worked out 
for one of the configurations should apply to the others. 

On the basis of the aforementioned comparison, if tube-portion 
weights are comparable, it would appear that other factors that are not 
related to heat transfer may ultimately fix the choice of the fin and 
tube configuration. Another utility of the present report relates to 
the presentation of a correct method of approach, both for black and for 
nonblack surfaces. 


Heat Transfer laboratory 
Department of Mechanical Engineering 
University of Minnesota 
Minneapolis, Minn. 

July 10, 1962 
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APPENDIX - ANALYSIS FOR NONBLACK SURFACES 

Consideration will be given here to the derivation of equations 
governing energy transport in the closed-sandwich configuration. The 
corresponding equations for the other configurations can be derived 
from those of the closed sandwich in a manner identical to that employed 
in the body of the report. Therefore, the other configurations will not 
be treated here. 

The radiant energy leaving a nonblack surface consists not only of 
emission, but also includes the reflected portion of the incident radi- 
ation. For diffuse reflection and emission at a surface, the reflected 
and emitted energies have the same directional distributions, and no 
distinction need be made between them. It is therefore convenient to 
define a new energy quantity, the radiosity B, as the sum of the emit- 
ted and reflected energies per unit time and area: 

B = eaT 4 + pH (Al) 

where H is the incident energy per unit time and area. For a gray 
surface, the reflectivity p is equal to (l - e), and additionally 
a = e. 


By utilizing the radiosity, the analysis previously given for the 
black- surfaced, closed-sandwich configuration can be readily generalized. 
It will be assumed here for simplicity that all surfaces are gray and 
have the same emissivity, but the latter restriction can be lifted with- 
out difficulty. The energy balance (eq. (2)) on an elementary volume 
of fin as shown in figure 2 continues to apply, as does the net conduc- 
tion given by equation (3). The net radiation is equal to the difference 
between the energy emitted and energy absorbed. The emitted energy is 

£ecrT 4 (x)J 2 dA x , while the absorbed radiation is aH(x)dA x . With this, 
the expression for the net radiation is 

(®W) ra4 - €[2 cT 4 (x) - < A2 > 

The incident energy H(x) comes from two distinct zones as before: 
(1) from the opposite fin and (2) from the base surfaces. The first of 
these contributions has essentially the same mathematical representation 
as before, except that the emissive power oT 4 (y) is replaced by the 
radiosity B(y). Thus, the generalization of equation (5) applicable to 
nonblack surfaces is 





B Cy) di ’dA x -dA y 


(A3) 
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The contribution of the base surfaces to H(x) must be considered in 
somewhat greater detail. Whereas an isothermal black surface has uni- 
form emissive power, it is not necessarily true that an isothermal non- 
black surface has uniform radiosity. In general, the radiosity on such 
a surface will not be uniform. Thus, for the case of nonblack surfaces, 
it is necessary to find the contribution to H(x) from each element of 
base surface and then integrate to find the total contribution. This is 
identical to the procedure employed in deriving expression (5). Pro- 
ceeding along the same lines, but replacing emissive powers by radio- 
sities, the contribution of the base surfaces to H(x) is found to be 





B ( z l) dF dA x -dA 1 



B ( z 2 ) c ^ B ’dA x; -dA2 


(A4) 


The radiosities B(z^) and B(zg) are equal when z^ = Z 2 because of 
symmetry. 

The incident energy H(x) is the sum of expressions (A3) and (A4). 
These may be employed to evaluate equation (A2) for the net radiation 
at dA^ With this and with the net conduction from equation (3), the 
heat balance (eq. (2)) becomes 


^ & = 2oT* - 
e dx 2 


y*2L 

«'y=0 

i 


dA^-dAy 


B(zi)dP d A x _dA 1 + 


h 

/ 

%/ z 9 =0 


B ( z 2) dB ’dA x -dA, 


(A5) 


It is convenient to introduce dimensionless variables X, Y, and Q as 
defined by equation (9), and in addition 


P = B/aT^ Z-l = z-j/h Z 2 = z 2 /h 


(A6) 
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Also, because of symmetry B(y) = B(2L - y), B(z^) = B(zg), B(z-j_) - 
B(h - Zj_). Incorporating all this into equation (A5), there is obtained 



in which the angle factors between dA^. and dA-^ and between dA x and 
dAg have been evaluated by application of equation (12) with dy re- 
placed by dz. The boundary conditions on 9 remain as stated in equa- 
tion (ll). 

Inspection of equation (A7) reveals that there are three unknown 
functions: the temperature distribution 9(x) and the radiosity dis- 

tributions p(Y) and p(Z). Two additional relations between the un- 
knowns are therefore required. These relations may be derived by making 
use of the radiosity equation (Al). First, applying this to the location 
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x and noting that the incident flux h(x)/ctT^ is given by the quantity 
within the braces on the right side of equation (A7), there follows 

P(X) = €0 4 + (1 - e) ^Terms within braces of eq. (A7 ) ) (A8) 

Because of symmetry, 0(x) = 0(Y) when X = Y. So, p(x) and 0(y) are 
the same functions; only the independent variables have been inter- 
changed. It then follows that equation (A8) is essentially an integral 
equation for 0(x). An additional relation between the unknowns is ob- 
tained by applying the radiosity equation to a surface element on one 
of the base surfaces. After evaluating the incident radiation from the 
fins and the opposite base surface, equation (Al) becomes 
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Because of symmetry, P(Z^) = PCZg) when Zj = Zg. So, 0(Zp and. P(Zg) 
are the same functions; only the independent variables have been inter- 
changed. Equation (A9) is therefore an integral equation for finding 

P(Zi)- 

The governing equations for the energy transfer in the closed- 
sandwich configuration with nonblack surfaces are embodied in (A7), 

(A8), and (A9)* Simultaneous solution of this nonlinear, integro- 
differential system is required to determine the heat-transfer charac- 
teristics, By comparing this mathematical system with equation (17), it 
is clear that the consideration of nonblack surfaces adds a computational 
complexity to the problem in addition to new unknowns and parameters* 
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Figure 5- - Heat loss due to conduction in fin relative to total heat loss. Closed- sandwich 
c onf igur at ion* 
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Figure 9. - Heat-loss results for central-fin configuration. 







Figure 11. - Concluded. Heat-loss comparison on the basis of equal fin weight. 
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Figure 12. - Comparison of results from superposition method and from numerical solutions 












